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ABSTRACT 

New far-infrared and sub-millimeter photometry from the Herschel Space 
Observatory is presented for 61 nearby galaxies from the Key Insights on 
Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. 
The spatially-integrated fluxes are largely consistent with expectations based on 
Spitzer far-infrared photometry and extrapolations to longer wavelengths using 
popular dust emission models. Dwarf irregular galaxies are notable exceptions, 
as already noted by other authors, as their 500 /im emission shows evidence for 
a sub-millimeter excess. In addition, the fraction of dust heating attributed to 
intense radiation fields associated with photo-dissociation regions is found to be 
(21 ± 4)% larger when Herschel data are included in the analysis. Dust masses 
obtained from the dust emission models of Draine Sz Li are found to be on average 
nearly a factor of two higher than those based on single-temperature modified 
blackbodies, as single blackbody curves do not capture the full range of dust 
temperatures inherent to any galaxy. The discrepancy is largest for galaxies 
exhibiting the coolest far-infrared colors. 

Subject headings: ISM: general — galaxies: ISM — infrared: ISM 



1. Introduction 

The amount of dust contained within a galaxy reflects the integrated chemical enrich- 
ment of the interstellar medium through accumulated episodes of star formation and any 
merger/accretion events, coupled with the galaxy's history of dust grain formation and de- 
struction. Thus, the total amount of dust observed within a galaxy can be tied to its star 
formation history. From an observational viewpoint, estimating a galaxy's dust mass de- 
pends critically on far-infrared/sub- millimetei0 photometry and the inferred distribution of 
dust grain temperatures; long wavelength data are crucial for probing cool 15-20 K dust with 
any accuracy, and dust emitting at thi s temperature range makes up the bulk of the dust 



mass in a typical star- forming galaxy (IDunne &: Ealesl 1200 ll ). Recent observational efforts 



have shown evidence for excess emission at sub-millimeter wavelengths in galaxies, particu- 
larly in low-metallicity systems such as the Small Magellanic Cloud and other dwarf galaxies 



Copenhagen, Denmark 

^'Far-infrared" and "sub-millimeter" are denned here as 40 < A < 300 and 300 < A < 1000 ^m, 
respectively. 
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(e.g.,|B o 
2011al lbl: 



atto et al. 2000 



Galametz et al. 



Galliano et al. 2005 



2011 



Bot et al. 2010: Gordon et al. 2010: Ade et al. 



Galliano et al 



201ll ). If the excess emission is interpreted as 



evidence for particularly cold dust (e.g., T < 10 K), substantial amounts must be present in 
order to produce the observed emission, and the resulting dust-to-gas mass ratios inferred 
from these observations are s ignificantly higher than what would be expected based on a 



galaxy's metallicity (see, e.g., iGalametz et al.l 120091 ; iMeixner et al.ll2010l : lAde et al.ll2011al ). 



Alternative explanations for the excess emission include resonances due to impurities in the 



size and/or comi 


position (see 


Lisenfeld et al. 


2002 




Ag 


uirre et al. 


Menv et al. 


2007 




Galametz et al. 


2009 




Meixner et al. 


2010 


, and 



references therein). 

Understanding the physical origin of the sub-millimeter emission in nearby galaxies obvi- 
ously requires sensitive sub-millimeter data. However, it has been difficult to obtain reliable 
sub-millimeter photometry. The opacity of the atmosphere limits ground-based observations 
to a few partially-transparent sub-millimeter spectral windows, and previous space-based 
missions such as COBE, IRAS, ISO, Spitzer, and AKARI lacked the sensitive instrumen- 
tation and the requisite angular resolution fo r resolved studies of nearby galaxies at these 
wavelengths. The Herschel Space Observatory (IPilbratt et al.ll2010f ) allows an unprecedented 
look into the long wavelength emission from galaxies. Herschel provides impressive angular 
resolution in the far-infrared/sub- millimeter wavelength range and unparalleled sensitivity 
to low surface brightness emission at 250, 350, and 500 /im, wavelengths that bridge an im- 
portant spectral gap between space-based galaxy surveys carried out at 160-170 /im (ISO, 
Spitzer, AKARI) and those from the ground at 850-870 /im; there are precious few de- 
tections of galaxies at 350 or 450 am of JCMT/SCUB A or CSO/SHARC heritage (e.g., 
Benford et al.lll999t iDunne fc Ealesll2001t iKovacsl 120061 ) . Filling in this gap in wavelength 
coverage is particularly important for the coldest galaxies in the local universe, in which the 
dust emission peaks at ~150 /im or longer wavelengths. 



In this effort we present Herschel PACS (jPoglitsch et al.ll2010l ) and SPIRE ( IGriffin et al. 



20101 ) broadband photometry for the 61 galaxies in KINGFISH (Key Insights on Nearby 
Galaxies: A Far- Infrared Survey with Herschel). The flux densities provided here are 
spatially-integrated or "global" values, and thus they do not directly take advantage of one of 
the key features of Herschel data: angular resolution. However, the sharpness of the imaging 
comes into play by enabling superior identification and removal of contamination from neigh- 
boring and background galaxies. While this study of global flux densities only indirectly 
profits from HerscheFs superior resolution, othe r studies inve s tigate the spatially-resolve d 
pro perties of KINGFISH galaxies. For example, IWalter et al.l ( 1201 ll ). lAniano et al.l ( 120121 ). 
and iGordon et al.l (j2012af ) explore how the infrared spectral energy distributions vary with 
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l ocatio n in KINGFISH targets, iGalametz et al.l (120121). Iffinz et all (120121). and Hunt et al. 



Crocker et al. 



(boi2h 



(120121 ) study cold dust emission in the outer disks, iLi et al.l (120121 ) and 
characterize the infrared properties of H II regions, etc. 

The global flux densities presented here are used to see how well the far-infrared/sub- 
millimeter emission i n nearby gala x ies ca n be modeled using the theoretical spectral energy 
distribution curves of iDraine fc Lil ( 120071 ). curves that are based on the dust emission prop- 
erties of the Milky Way. We also quantify whether the sub-millimeter data provide new 
insight into galaxy spectral energy distributions or whether their values are expected based 
on extrapolations from shorter wavelength data. Finally, we search for evidence for a global 
sub-millimeter excess in KINGFISH galaxies, explore the unique characteristic s of the low- 
metallicity systems in the sample, and contrast the dust masses found through IDraine &: Li 
( 120071 ) fits with those from the common approach of fitting single-temperature modified 
blackbodies. 

Section [2] describes the sample and § [3] reviews the observations as well as the data 
processing and aperture photometry procedures that have been adopted here for KINGFISH 
imaging. In § 0] the spatially-integrated photometry data are presented in addition to the 
results of fits to the observed spectral energy distributions, and § [5] provides a summary of 
our findings. 



2. Sample 



The KINGFIS H sample of galaxies draws from the Spitzer Infrared Nearby Galaxies 
Survey (SINGS; see iKennicutt et al.ll2003l ); fifty-seven SINGS galaxies are in the KINGFISH 
sample, plus NGC 2146, NGC 3077, NGC 5457, and IC 342. The 61 galax ies were selected 



to sp an wide ranges in luminosity, optical/infrared ratio, and morphology (IKennicutt et al. 
20121 ). The luminosity ranges over a factor of 10 4 (but all are sub-LIRG, or nearly so in 



the case of NGC 2146, with L IR < 10 L ), the range of optical/infrared ratios covers a 
factor of 10 3 , and all "normal" galaxy types are represented. There are several galaxies with 
nuclei that are clearly distinguished by Seyfert or LINER characte ristics, but none of th e 
galaxies has their global luminosity dominated by an active nucleus ( jMoustakas et al.ll2010l ). 
The sample also spans representative ranges in metallicity, gas fraction, H I/H 2 ratio, star 
formation rate, and bar strength. 



3. Observations and Data Processing 



The ob servational program and d ata processing procedure s for KINGFISH are describe d 
in detail in lEngelbracht et al.l ( 120101 ) , ISandstrom et al.l ( 120101 ) , and iKennicutt et al.l ( 120121 ) . 
A brief description is provided here. All 61 KINGFISH galaxies were imaged with PACS 
and SPIRE. 



3.1. PACS Data 



PACS imaging was obtained in scan mode, along two perpendicular axes for improved 
image reconstruction, at the medium scan speed of 20" s _1 . The 45° orientation of the 
array with respect to the scan direction contributes to a more uniform spatial coverage. 
Two Astronomical Observation Requests (AORs) were carried out for joint 70 and 160 /im 
imaging, and an additional two AORs were utilized for joint 100 and 160 /im observations, 
resulting in a total of four AORs for 160 /im imaging. Three or six repetitions were carried out 
for each AOR, depending on an individual galaxy's far-infrared surface brightness as gauged 
from Spitzer/MIPS data. The integrations achieved per pixel lead to approximate la surface 
brightness sensitivities of cr s k y ~ 5, 5, and 2 MJy sr _1 at 70, 100, and 160 /im, respectively, for 
the fainter subset of galaxies and approximately \/2 times larger for the brighter subset. The 
PACS calibration uncertainties are e ca \ >v / f v ~ 5%, according to Version 2.3 (08 June 2011) 
of the PACS Observer's Manual. 



The raw ("Level 0") data were processed using Version 5.0 of HIPE (IOttll2010l ). Be- 
sides the standard pipeline procedures, the conversion from Level to Level 1 data in- 
cluded second-le vel deglitching and corrections for any offsets in the detector sub-matrices. 
Scanamorpho3§ (Roussel 2011 ) Version 12.5 was used to process the Level 1 PACS scan map 
data. Its main task is to subtract the brightness drifts caused by the low-frequency noise 
(comprising both the thermal drifts of the telescope and detectors and the flicker noise of 
the individual bolometers), before projecting the data onto a changeable spatial grid. The 
algorithm employs minimal assumptions about the noise and the signal, and extracts the 
drifts from the data themselves, taking advantage of the redundancy built in to the scan 
observations. With the nominal settings used by the KINGFISH survey, the drifts can be 
determined on timescales greater than or equal to 0.7 s at 70 and 100 /im, and 0.9 s at 
160 /im (for a sampling interval of 0.1 s). These timescales correspond to lengths between 
1.5 and 2.5 times the beam FWHM, from 160 /im to 70 /im. Second-level deglitching was 



2 http: / /www2. iap.fr / users / roussel/herschel 
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performed, and the option to detect and mask brightness discontinuities was also used. The 
data are weighted by the inverse square high-frequency noise of each bolometer in each scan. 

The ( "Level 2" ) output of Scanamorphos is in the form of a FITS cube for each filter. 
The four planes are the signal map, the error map, the map of the drifts that have been 
subtracted, and the weight map. There is currently no propagation of the errors associated 
with the successive processing steps in the pipeline. In each pixel the error is defined as the 
unbiased statistical estimate of the error on the mean. The brightness unit is Jy pixel -1 , 
and the pixel size is ~ one-fourth of the beam FWHM, i.e., 1.4" at 70 /im, 1.7" at 100 /im, 
and 2.85" at 160 /im. 



3.2. SPIRE Data 

SPIRE imaging data were taken in Large-Map mode, to an extent tailored to each 
galaxy's size (out to at least ~ 1.5 times the optical size). Either two or four scans were 
obtained for each galaxy, based on its Spitzer/MIPS far-infrared surface brightness. The re- 
sulting la limiting surface brightnesses are approximately a s k y ~ 0.7, 0.4, and 0.2 MJy sr _1 
at 250, 350, and 500 /im, respectively, for the fainter subset and \/2 larger values for the 
brighter galaxies. Calibration uncertainties for SPIRE data are estimated at e c& \/ f v ~ 7%, 
following Version 2.4 (07 June 2011) of the SPIRE Observer's Manual. However, the uncer- 
tainties are strongly correlated between the three bands and thus the uncertainty on some 
SPIRE quantities such as the f v (250fim)/f u (500fim) flux density ratio, for example, are less 
than the simplistic \/2 x 7% expectation. 

SPI RE observations fo r six of our galaxies were obtained in the Herschel Reference 



Survey flBoselli et all 120101 ): NGC 4254, NGC 4321, NGC 4536, NGC 4569, NGC 4579, 
and NGC 4725; those observations were not duplicated for KINGFISH. The only notable 
difference between the SPIRE observations for the Herschel Reference Survey and those for 
KINGFISH is that three scans were employed (versus either two or four for KINGFISH 
observations, as described above). 

The raw SPIRE data are processed through the early stages of HIPE (Version 5.0) to 
fit slopes to the data ramps and to calibrate the data in physical units. A line is fit to 
the data for each scan leg after masking out the galaxy, and this fit is subtracted from 
the data. Discrepant data (usually due to a rogue bolometer, of which there are < 1 per 
map) are also masked, and the data are mosaicked using the mapper in HIPE. The map 
coordinates are then adjusted so that the position of the point sources (measured using 



StarFinder; iDiolaiti et al.ll2000f ) match those in the MIPS 24 /im images. Finally, the images 
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are converted to surface brightness units by dividing by the beam areas published in the 
SPIRE Observer's Manual: 423, 751, and 1587 □" at 250, 350, and 500 /mi, respectively. 
Pixel sizes are 6", 10", and 14" at 250, 350, and 500 /mi, respectively. 



3.3. Sky Subtraction 

At far-infrared/sub- millimeter wavelengths the emission from the sky (above the at- 
mosphere) largely comes from Milky Way cirrus and background galaxies. However, the 
bolometer arrays of Herschel are not absolute photometers, and thus any map produced by 
any software is the superposition of an estimate of the true sky emission and an unknown 
(large) offset. Hereafter this superposition is referred to as simply the "sky". While the 
post-pipeline processed KINGFISH SPIRE and PACS images have their overall sky levels 
removed to zeroth order, a procedure has been adopted to remove a more refined local sky 
value for each galaxy. To accomplish this local sky subtraction, for each PACS and SPIRE 
image a set of sky apertures has been defined that collectively circumscribe the galaxy, pro- 
jected on the sky close enough to the galaxy to measure the "local" sky but far enough 
away to avoid containing any galaxy emission (Figured]). The emission from any prominent 
neighboring and/or background galaxies that are projected to lie within the sky apertures 
is removed before the sky is estimated. The total sky area, derived from the sum of the ar- 
eas from all sky apertures, is typically significantly greater than that covered by the galaxy 
aperture itself, thereby limiting the contribution of uncertainty in the sky level to the overall 
error budget. The mean sky level per pixel is computed from the collection of these sky 
apertures, the value is scaled to the number of pixels in the galaxy aperture, and the result 
is subtracted off from the overall galaxy aperture counts (all done within IRAF/IMCNTS). 



Aniano et al.l (120121 ) follow a different procedure for subtracting the sky emission from 



KINGFISH imaging, including fitting a tilted plane t o the s ky for each galaxy instead of 



a single value approach adopted here. I Aniano et al.l (120121 ) study the spatial variations 



in the far-infrared/sub- millimeter spectral energy distributions and thus a more detailed 
characterization of the local sky is necessary. The effects of most sky gradients cancel out in 
extracting spatially integrated fluxes; the two approaches yield generally consistent global 
fluxes. 
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3.4. Aperture Photometry 

The elliptical apertures used for global photometry are listed in Table [TJ The apertures 
are chosen by eye to encompass essentially all of the emission at every wavelength; aperture 
corrections described below are incorporated to recover the amount of any light that lies 
beyond these apertures. The average ratio of aperture major axis length 2a to the de 
Vaucouleurs D25 optical major axis is 1.45 (with a la dispersion in this ratio of 0.45). The 
same aperture is used to extract the flux at each wavelength st udied here, and t hey are very 



similar, if not identical, to those used for SINGS photometry (IDale et al.l 120071 ). for the 57 
galaxies that overlap the two samples. As a test of the robustness of the aperture choices, 
the global flux densities using these apertures are compared to the values obtained using 
apertures with 5% larger and 5% smaller semi-minor and semi-major axes. The impact of 
using ~10% larger or ~10% smaller aperture areas is a median difference of less than 1% on 
the flux densities for all wavelengths. 

Prior to extracting fluxes from aperture photometry, any emission from neighboring 
or background galaxies is identified and removed from the area covered by each aperture. 
The identification is assisted by ancillary data at shorter wavelengths and higher spatial 
resolution (e.g., Spitzer/IRAC 3.6 and 8.0 /zm, HST optical, and ground-based Ha imaging). 
The removal is accomplished via IRAF/IMEDIT by replacing the values of contaminated 
pixels with the values from a random selection of nearby sky pixels, thereby incorporating 
the same noise statistics as the sky. Usually the removal of such emission affects the global 
flux at less than the 1% level, but in a few cases the impact is quite important, e.g., NGC 1317 
lies within the aperture of NGC 1316, and background galaxies in the fields of the fainter 
dwarfs like Ho II and DDO 053 would c ontribute significan tly to the integrated flux (by up 



to ~ 30 — 50%) if not removed (see also I Walter et al.l 120071 ) . 



Diffraction inevitably results in a small portion of the galaxy emission appearing beyond 
the chosen apertures, however, and thus aperture corrections are formulated to mitigate this 
effect. Aperture corrections are empirically determined from a comparison of fluxes from 
smoothed and unsmoothed Spitzer/IRAC 3.6 /zm imaging, which has a native resolution of 
~l'/5. The aperture correction for a given PACS or SPIRE flux is the ratio of the flux from 
the unsmoothed 3.6 /j,m image to the flux from the 3.6 /xm image smoothed to the same PSF 
as the Herschel image in question. Due to the typically generous aperture size and sharp 
Herschel PACS and SPIRE PSFs, the amplitudes of the KINGFISH global photometry 
aperture corrections are typically quite small, with median values of 1.0 at all wavelengths 
and maximum values of 1.03 for PACS and between 1.07 and 1.13 for SPIRE. This technique 
assumes that a galaxy's profile in the far-infrared matches that of its (mostly) stellar profile 
in the near-infrared, and there may in fact be appreciable differences in the two emission 
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profiles. 

The uncertainties in the integrated photometry e to tai are computed as a combination in 
quadrature of the calibration uncertainty e ca i and the measurement uncertainty e s k y based 
on the measured sky fluctuations and the areas covered by the galaxy and the sum of the 
sky apertures, i.e., 

6 total = ^cal + C sky (!) 

with 

Csky = Csky^pix \J N pix + N pix 2 / N sky (2) 

where <7 s k y is the standard deviation of the sky surface brightness fluctuations, f2 P i x is the 
solid angle subtended per pixel, and iVp ix and A^ sky are the number of pixels in the galaxy 
and (the sum of) the sky apertures, respectively. For the few sources undetected by Herschel 
imaging, 5a upper limits are derived assuming a galaxy spans all N pix pixels in the aperture, 

fv(5a upper limit) = 5 e s k y . (3) 



Results 



4.1. Flux Densities 



Table [2] presents the spatially- integrated flux densities for all 61 KINGFISH galaxies for 
all six Herschel photometric bands. The tabulated flux densities inc 



ude aperture corrections 

(§ 13.41) and have Galactic extinction (ISchlegel. Finkbeiner. &: Davislll998h removed assuming 
AyjEiB — V) ~ 3.1 and the reddening curve of iLi fc Draine! feoOlh l 3 ! No color corrections 
have been applied to the data in Table [2j The most recent calibrations are used for both 
SPIRE and PACS photometry, including the "FM,6" PACS calibration which lowers the 
fluxes for extended sources by 10-20% compared to the previous calibration, as described in 
the HIPE 7.0.0 documentation. 

The superior sensitivity and angular resolution of Hersche l enables a more detailed 
investigation of the faintest galaxies in our sample. For example, iDale et al.l (120071 ) provide 
marginally significant MIPS flux densities at 70 and 160 fim for NGC 1404 and DDO 165, but 
they caution that the emission appearing within the apertures for these galaxies potentially 



3 The corrections for Galactic extinction are very small at these wavelengths, with the largest correction 
being 0.4% for IC 0342 at 70 /im, which lies at a Galactic latitude of +10° behind a foreground veil of 
E(B - V) ~ 0.56 mag. 
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derives from background galaxies. It is now clear based on the Herschel maps that these 
targets were indeed not detected by Spitzer (nor Herschel) at A > 70 /im. 

The suite of far-infrared filter bandpasses available for Herschel and archival Spitzer 
data allows a direct comparison of the global flux densities measured for the SINGS and 
KINGFISH galaxy samples. The Spitzer and Herschel fluxes at 70 and 160 /im are on average 
consistent; Figure [2] shows the (error-weighted) ratios of Spitzer/ Herschel flux densities agree 
fairly well, after accounting for differences in the Spitzer and Herschel spectral responses 
and calibration schemes. The agreement at 70 /im is within 3% (with a 12% dispersion in 
the ratio), while at 160 /im the MIPS flux densities are typically 6% larger (with a 16% 
dispersion). Galaxies fainter than ~1 Jy show a much larger dispersion in these ratios, but 
the flux densities for these targets are more susceptible to errors in sky estimation (§ 13. 3p . 

Figures [3] and H] provide a color-color snapshot of the Herschel global photometry. As ex- 
pected, a clear correlation is seen in Figure [3] when the flux density ratios on both axes involve 
wavelengths that straddle the broad infrared peak of emission for most galaxies. Figure HJ 
on the other hand, demonstrates that the galaxy spectral energy distributions in general do 
not form a simple one-parameter sequence. Galaxies are more complicated, with mixtures 
of dust temperatures and distributions of grain properties that vary from one galaxy to an- 
other. The galaxy types are fairly well distributed in terms of their infrared/sub- millimeter 
colors, though the Sc and Sd spirals tend to cluster toward cooler infrared colors (i.e., smaller 
values of / I/ (70/im)// iy (160/im), / I/ (70/xm)// iy (250/xm), and / I/ (100/im)// 1/ (500/xm)) and the 
Magellanic Irregulars (type Im) have relatively large 500 /im flux densities. This latter issue 
will be revisited in § 14.61 



4.2. The Observed Spectral Energy Distributions 



Figure [5] shows the observed infrared/sub- millimeter spectral energy distributions for 
the KINGFISH sample. Included in each panel, when available, are the 2MASS K s , ISO 
6.75 and 15 /mi, Spitzer 3.6, 4.5, 5.8, 8.0, 24, 70, and 160 /im, IRAS 12, 25, 60, and 
100 /im, Herschel 70, 100, 160, 250 350, and 5 00 /im and SCUBA 450 and 850 /im band 
fluxes, derived from this work and iDale et al.l (120071 . 120091 ) . These data nomina lly reflect 
the global emission at each wavelength, but as pointed out in iDraine et al.l (120071 ). a subset 
of the SCUBA images suffers from various technical and observational issues. The data 
processing for scan mapped SCUBA observations (NGC 4254, NGC 4579, and NGC 6946) 
removes an unknown contribution from extended emission; the areas mapped by SCUBA for 
NGC 1097, NGC 4321, an d NGC 4736 were small and thus any errors in the large aperture 
corrections determined by IDale et al.l (120071 ) for these three systems would have a significant 
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impact on their inferred global fluxes; and the extra-nuclear sub-millimeter emission at 
850 /im is unreliably mapped in NGC 4594 due to contamination by an AGN. Indeed, for 
all these special cases except NGC 4736, the SCUBA data appear to fall appreciably below 
expectations based on extrapolations from the superior Herschel data. 



4.3. Fits to the Observed Spectral Energy Distributions 



To extract physical parameters from th e broadband spectr al data, the spectral energy 
distributions were fitted with the models of iDraine &: Lil (120071 ). models based on mixtures 
of amorphous silicate and graphitic dust grains that effectively reproduce the average Milky 
Way extinction curve and are consistent with observat i ons o f PAH features and the vari- 
ety of infrared continu a in lo cal galaxies. IDraine fc Lil (120071 ) use the size distributions of 
Weingartner fc Draind (120011 ) for dust in the diffuse Milky Way, except for ad j ustm ent of 
the parameters that characterize the PAH size distribution. The IDraine &: Lil (120071 ) dust 
mod els use the far - infrar e d and sub-mi ll imete r opacities for graphite and amo rphous silicate 
fromlLi fc Draind (120011 ). iLi fc Draind (120011 ) used the graphite opacity from IDraine &: Lee 
(119841 ). but made small modifications to the amorphous silicate opacity. The imaginary 
part of the amorphous silicate dielectric function e 2 (A) was adjusted in order for the model 
to better matc h the average high Galactic latitude dust emission spectrum measured by 
COBE-FIRAS dWright et al.lll99ll : [Reach et ali fl~9~95l iFinkbeiner et al.l fl~9~9~9l ) . The adjust- 
ments were modest: 62(A) was unchanged for A < 250 /im, and modified by less than 12% 
for 250 / mi < A < 1100 /xm. With this dielectric function for the amorphous silicate compo- 



nent, the IDraine fc Lil (120071 ) model gives generally good agreement with the o bserved sub 
millim eter emission from the Milky Way diffuse interstellar medium. Thus the IDraine fc Li 
(120071 ) model has in effect been "tuned" to reproduce the diffuse emission fr om the local 
Milky Way. While the dust model used here is referenced as coming from IDraine &: Li 
(120071 ). in fact two small changes have been incorporated since that publication: i) there 
have been some small changes in some of the PAH band parameters, and ii) the graphite 
dielectric function has b een modified t o bro aden out an opacity peak near 30 /im. These 
changes are described in lAniano et al.l (|2012| ). 



Building upon an idea put forth by iDale et al.l ( 120011 ). IDraine fc Lil ( 120071 ) model inter- 
stellar dust heating within a galaxy with a 5-function in interstellar radiation field intensity 



U coupled with a power-law distribution [/ min < U < U max , 



dM dust /dU = M dust 



(l- 7 )5(t/-t/min)+7-T 



a — I 



Ut 



JJl-a 



-ET 



(4) 



where U is normalized to the local Galactic interstellar radiation field, dMA mf is the differ- 
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ential dust mass heated by a range of starlight intensities [U, U + dU], Md us t is the total dust 
mass, and (1 — 7) is the portion of the dust heated by the diffuse interstellar radiation field 
defined by U — U m [ n . The minimum and maximum interst ellar radiat i on fie ld intensities 
span 0.01 < U min < 30 and 3 < log£/ max < 8. See § 5.5 of iDale et al.l (120011) for a ph ysi- 



cal motivation of the power law distribution in U, and Figure 4 of iDraine fc Lil (120071 ) for 
examples of translating U to dust temperature for different grain sizes. 

A sum of three different spectr al energy dist r ibutio ns is fit to each galaxy: a blackbody 
of temperature T* = 5000 K, which ISmith et al.l (120071 ) find to be a good appr oximation to 



the st ellar profile beyond 5 /im, along with two dust components. Following IDraine et al. 
(120071 ). the sum can be expressed as 



model 



^;f; -By f-^ - * 



M, 



dust 



4nD 2 



PAH; f^min) !Pu(qpAu, U min , U max , a)] , (5) 



where fl* is the solid angle subtended by stellar photospheres, D is the distance to the galaxy 
and 7 and (1 — 7) are the fractions of the dust mass heated by the "power-law" and "delta- 
function" starlight distributions, respectively, pi ((/pah, Umin) and Pv(lpAa, U^, U masL , at) 
are, respectively, the emitted power per unit frequency per unit dust mass for dust heated 
by a single starlight intensity £7 mm , and dust heated by a power-law distribution of starlight 
intensities dM/dU oc U~ a extending from U m in to U max . The U = U m i n component may be 
interpreted as the dust in the general diffuse interstellar medium. The power-law starlight 
distribution allows for dust heated by more intense starlight, such as in the intense pho- 
todissociation regions (PDRs) in star-forming regions. For simplicity, emission from dust 
heated by U > U m i n will be referred to as the "PDR" component, and the emission from 
dust heated by U — U min will be referred to as the "diffuse ISM" component. Finally, the 
fractional contribution by total dust mass from PAHs, denoted gpAH, varies between 0% and 
12% with a model grid spacing of 0.1% in gpAH- 



Draine et al.l (120071 ) find that fits to the (SINGS) global spectral energy distributions 



of nearby galaxies are insensitive to the minimum radiation field intensity, the maximum 
radiation field intensity, and the power-law parameter a. We a dopt their choice t o fix 
U max = 10 6 and a = 2 to minimize the number of free parameters. IDraine et al.l (120071 ) use 
a minimum value of 0.7 for U m i n , but we choose to extend this range down to 0.01 due to 
the availability of SPIRE data longward of 160 /mi and the resulting potential for having 
detected very cold dust emission. The free parameters fl*, M dust , qpAR, f^mm, and 7 are found 
via x 2 minimization: 



X 



E 



/ -fobs 



/ model \ 2 
Jv,b I 



a 



obsN2 



model \ 2 ' 



(6) 



where /™ 6 ° is the model flux density obtained after convolving the model with the b filter 
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bandpass, cr£ bs is the uncertainty in the observed flux density, and <j™ odel is set to 0.1/™ fe odel 
to allow for the uncertainty intrinsic to the model. 



Figure displays the fits of the iDraine &: Lil ( 120071 ) models to the combined broadband 
observations from the Spitzer and Herschel observatories. The median reduced chi-squared 
value is near unity (~ 0.7), and with just a few exceptions the fits are quite reasonable. The 
most challenging spectral energy distributions to fit involve spatially variable Milky Way 
cirrus coupled with a faint target, and thus any errors in determining the value of the local 
sky has a relatively large impact on the inferred fluxes (i.e., dwarf galaxies such as DDO 053, 
M81 dwarf B, and the faint elliptical NGC 0584). 

A wealth of information can be extracted from such fits. Figure for example, uses 
these fits to provide a glimpse into how the global spectral energy distributions depend on the 
star formation rate and total infrared luminosity. The infrared spectral energy distributions 
typically peak at shorter wavelengths for KINGFISH galaxies with higher star formation 
rates and infrared luminosities. There are exceptions to these generalizations, however, 
especially for lower luminosity systems. A full tabula t ion of the output parameters for the 
KINGFISH sample will be presented in lAniano et al.l (120121 ) . Here we restrict our analysis 
of the output parameters to i) evaluating the impact of i ncluding the Hersc hel data in these 
fits, and ii) comparing the dust masses found through IDraine &: Lil (120071 ) fits with those 
from single-temperature modified blackbody fits. 



4.4. Spectral Energy Distribution Fit Parameters 

Figure [7] compares (ratios of) the output parameters 7, ^pah, f^min, and M dust when the 
fits are executed with and without the inclusion of Herschel photometry. All four parameters 
are relatively unchanged, on average, when Herschel broadband data are added to those from 
Spitzer. The largest average deviation in the ratio from unity is seen in the top panel, where 
the fraction of dust heated by photo-dissociation regions found by using both Spitzer and 
Herschel data is on average (21 ±4)% larger than that using just Spitzer data. Interestingly, 
the largest dispersions in the distributions in Figure [7] are for U m i n and Md us t, indicating the 
importance of Herschel data in assessing these parameters. In addition, all four parameters 
show ratio distributions that are fairly evenly distributed about their means, though the 
scatter shrinks for cooler galaxies. At first blush it may be surprising that the inclusion 
of Herschel far-infrared/sub- millimeter has any impact on a parameter such as (/pah that 
is sensitive to mid-infrared PAH features, but recall that (/pah is the PAH mass abundance 
with respect to the total dust mass, and clearly Herschel photometry has an important role 
in determining the latter. Finally, even though U m { n was allowed to go as low as 0.01, the 
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smallest fitted value usi ng the combined Sy itzer and Herschel datasets is 0.6, similar to the 
U min floor advocated by iDraine et al.l ( 120071 ) . 



Figure E]shows the same ratios plotted as a function of oxygen abundance (IMoustakas et al. 



2010l ). T he dependence o f KING FISH dust masses on metallicity are consistent with those 
found by iGalametz et al.l (120111 ) in a study of 52 galaxies with sub-millimeter data: dust 
masses computed for metal- rich (metal-poor) galax ies are smaller (larger) when sub- millimeter 
data are included in the fit. IGalametz et al.l (120111 ) argue that most metal-rich galaxies have 
their dust emission peak in the far-infrared beyond 160 //m, and that sub-millimeter data 
are required to fine-tune dust measures for such systems. KINGF ISH galaxies with o xygen 
abundances 12 + log(0/if) greater (less) than 8.1 on the empirical iPilyugin &: Thuanl (120051 ) 
metallicity scale are computed to have an average of 0.06 ± 0.03 dex less (0.28 ± 0.09 dex 
more) dust mass when sub-millimeter data are used. This demarcation in metal abundance 
is similar to that studied by othe rs in quantifying, for example, the relative importance 
of PAH emission i n gal axies (e.g., iHunt et al.ll2005t IDraine et al.l 120071 ; ISmith et al.l 120071 : 
Engelbracht et al.l 120081 ). In short, perha ps KINGFISH data s how metallicity-dependent 



dust mass trends similar to those found by 
have more data to confirm any such trend. 



Galametz et al.l (120111 ). but it would be useful to 



4.5. Comparison with Dust Masses from Blackbody Fits 



Galaxy dust masses are typically estimated by fitting a single modified blackbody to a 
selection of infrared/sub- millimeter continuum fluxes, 



M, 



fuD 2 (i^Y 1 
~ k{vo)B v {T x ) \v) ' 



(7) 



or in some cases a superposition of two modified blackbodies of two different dust tempera- 
tures, 



M, 



dust 



fuD 2 
k(u ) 



+ (1-^(T 2 )(- 



(8) 



where k{pq) is the dust mass absorption coefficient at the reference frequency uq, T\ and T 2 
are the modeled dust temperatures, (3i and (3 2 are the dust emis sivity indexes, and < x < 1; 
some authors choose to fi x the dust emissivity index(es) (e.g., iDunne fc Ealesll200ll ; iKovacs 
20061 ; iPascale et al.ll2009l ) . While such approaches provide quick and simple routes to gauging 
the dust mass, they do not capture the full range of dust temperatures inherent to any galaxy. 
However, due to their popularity it is instructive to compare blackbody-based dust masses 
to those determined from more nuanced models. 
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Figure [9] compares the dust masses obtained by using a single-temperature modified 
bl ackbody (Equation !?!) to those obtained in Section 14.41 Both approaches utilize the results 
of iLi fc Draind (120011 ) for dust absorption cross sections, and in particular u = c/250/iin = 
1.20 THz and k{vq) ~ 0.48 m 2 kg -1 are used for the SPIRE 250 /xm band in determining 
the (modified) blackbody dust mass. In addition, both the dust temperature Td and the 
dust opacity coefficient j3 are allowed to freely vary in the fit for each galaxy; the ranges for 
the fitted values are quite reasonable given that KINGFISH does not contain a n y extr eme 
objects: 18 < T d < 40 K and 1.2 < p < 1.9 (Figure QUI see Iskibba et all fl201lh for 
similar results based on modified blackbody fits to KINGFISH targets). In order to avoid 
contributions from stochastically heated dust grains in the computation of the blackbody- 
based dust masses, the top panel of Figure [9] shows results when only Herschel photometric 
bands from 100 /im through 500 fim are included in the fits. Results are not significantly 
different when 70 jum data are included (bottom panel); the median ratio in the top (bottom) 
panel is 0.53 (0.46). Figure [9] indicates that single temperat ure (modified ) black body dust 
masses typically underestimate the values obtained through a iDraine fc Lil (120071 ) formalism 
by nearly a factor of two (~ 1.9), and there is a trend toward larger underestimates for 
galaxies exhibiting cooler far-infrared colors. 

Similar results are obtained after fixing /3 to either 1.5 or 2.0, except for the situation 
where the blackbody fits are carried out over the wider 70-500 /im wavelength baseline for 
= 2.0. In that case the fitted dust temperatures are lower in order to compensate for an 
overly steep emissivity dependence on wavelength, resulting in larger quantitie s of dust and 
only a 25% underestimate in t he dust mass compar e d to those obtai n ed fro m IDraine fc Li 
( 120071 ) , echoing the findings in iMagrini et al.l (120111 ) . iDunne fc Ealesl (120011 ) likewise find a 
factor of two deficiency for single blackbody-based dust masse s, in their case comp ared to the 
dust mass derived from two (modified) blackbodies (see also lSkibba et al.ll201ll ). Figure ITT1 
shows a primary reason for the discrepancy: even when limited to A > 100/im photometry, 
single-temperature blackbody fits overestimate the dust temperature, thus underestimating 
the dust mass. The single-temperature model does not account for the contribution of warm 
dust emitting at shorter wavelengths and the temperatures are driven towards higher values 
in the attempt to fit both the short and long wavelength far-infrared emission. This effect 
is accentuated for galaxies with cooler large dust grains whose emission peaks at longer 
infrared wavelengths. A more comprehensive and detailed comparison of various dust mass 
indicators is studied in iGordon et al.l ( ]2012bl ) . 
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4.6. Sub-millimeter Excess 



As described in § (TJ several studies of dwarf galaxies show significant excess emission at 
sub-millimeter wavelengths. Inspection of Figure [5] shows that only a few spiral and elliptical 
KINGFISH galaxies have 500 /im fluxes that are noticeably above the fitted model curves 
(e.g., NGC 3049 and NGC 5474). However, it is interesting that dwarf /irregular/Magellanic 
galaxies preferentially show this excess. There are a total of twelve KINGFISH galaxies of 
Type Im (Magellanic irregular), Type 10 (non-Magellanic irregular) or Type Sm (Magellanic 
spiral), a nd of the nin e with detections at 500 /im, eight show hints of 500 /im emission 
above the braine k Lfl (|2007h model fit (IC 2574, Holmberg I, Holmberg II, M81 dwarf B, 
NGC 2915, NGC 4236, NGC 5398, and NGC 5408; see also Figure EJ. Quantitatively, the 
observed 500 /im excess can be defined with respect to the model prediction at 500 /im, 
namely 

z//^(500/im) observed - ^/ ! ,(500 / um) modcl 



f(500//m) 



(9) 



i// v (500//in) mode i 

A dozen KINGFISH galaxies show £(500/xm) > 0.6, including all eight of the dwarf /irregular /Magellanic 
galaxies listed above. However, it should be noted that interpreting £(500/im) for these sys- 
tems is complicated by t he fact that they are typically faint in the far-infrared/sub- millimeter 
(e.g.. IWalter et al.ll2007l ) and thus their measured flux values are the least reliable. Nonethe- 
less, it is interesting that the lowest metallicity objects in the KINGFISH sample are the 
sources most likely to show a sub-millimeter excess and thus potentially harbor the coldest 
dust or have peculiar dust grain characteristics. A detailed analysis of the sub-millimeter 



excess in KI NGFISH galaxies i s being carried out by iGalametz et al.l (]201l[ ). iGordon et al. 
<j2012ah . and lHunt et aD ( j2012h . 



5. Discussion and Summary 



Spatially-integrated far-infrared and sub-millimeter flux densities from the Herschel 
Space Observatory are provided for the 61 objects in the KINGFISH sample of nearby 
galaxies. All but three galaxies are detected in the far-infrared by PACS and all but four 
galaxies are detected in the sub-millimeter by SPIRE. The (color-corrected) Herschel PACS 
70 /im global flux densities agree with Spitzer MIPS 70 /im counterparts to within 3% (with 
a 12% dispersion), on average; the MIPS 160 /im flux densities are typically 6% larger than 
the PACS 160 /im flux densities (with a 16% dispersion). 



Draine & Li 


(2007 


) and 


Draine et al. 


(2007) 



fit to the combined Spitzer and Herschel 3.6-500 /im dataset. The fits provide constraints 
on the total dust mass Md us t, the PAH mass fraction (/pah, and the characteristics of the 
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radiation fields that heat the dust including the fraction 7 of the dust mass that is located in 
regions with U > U m { n , and the complementary fraction 1 — 7 that is located in the general 
diffuse interstellar medium. A full tabulation of the fit parameters will be presented in 
Aniano et al.l (120121 ); analysis of the fit results here is restricted to comparisons between fits 
with and without inclusion of the Herschel data. In general, the fits to Spitzer+ Herschel data 
produce parameter values that are consistent, to within a factor of 2, with those when just 
Spitzer data are fitted. However, the KINGFISH galaxies with oxygen abundances less than 
12 + log(0/H) < 8 tend to show larger dust masses and smaller PAH mass fractions when 
SPIRE data are included in the fits. A similar characteristi c oxygen abundance has been 
noted in other studies of the PAH abundance in g alaxies (e.g.. lHunt et al.ll2005l ; iDraine et al. 



20071 ; ISmith et al.l 120071 ; lEngelbracht et al.ll2008l ). In addition, the fraction of the dust mass 
located in regions with U > t/ min , 7, is (21 ±4)% larger when Herschel data are included 
in the fits. For a = 2 the fraction of the total dust lum inosity contributed by regions with 
U > 100 is given by Equation 29 of braine feU (bo07h : 

7 in(rj max /io 2 ) 



f(L dmt ,U >10 2 



(10) 



(1 -7)(1 - U min /U mSLX ) +-f\n(U max /U min )' 
This parameter is (16 ±4)% larger when Herschel data are included in the fits. These subtle 
differences in the fits presumably reflect the unprecedented ability of Herschel to properly 
account for contributions from cold dust grains, grains that sustain their meager thermal 
emission through heating by the diffuse radiation field that permeates a galaxy's interstellar 
medium. 

The presence of an excess of emission in the sub-millimeter has been noted in the litera- 
ture, particularly for low-metallicity galaxies. Most KINGFISH galaxies are well modeled by 
spectral energy distributions consistent with emission curves from the Milky Way and nearby 
galaxies, without needing to invoke an additional cold dust component. However, eight of 
the nine dwarf /irregular /Magellanic galaxies with detections at 500 /j,m show evidence for 
significant e xcess of emission at this wavelength, at least with respect to the expectations 
based on the lDraine fc Lil ( 120071 ) model fits. These excesses, in fact, are the reason their dust 
masses are larger when Herschel data are included in the fits described above, assuming these 
excesses are attributable to increased quantities of very cold dust. It is unclear why low- 
metallicity dwarf irregular galaxies exhibit a propensity for conspicuous cold dust emission. 
In fact, their spectral energy distributions do not typically peak at longer wavelengths than 
is seen for the more metal rich galaxies; the KINGFISH dwarf galaxies are not colder than 
average, they simply show 500 /iin excesses. Perhaps such environments pr omote unusual 
dust emissivities that lead to the observed excesses (see also lBot et al.l I2010L for additional 
explanations). 



It is commonplace to find in the literature dust masses based on fits to single modified 
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blackbody profiles, with the dust temperature and dust emissivity modifier serving as 
potential free parameters. Blackbody-bas ed dust mas s es are on average a factor of ~ 1.9 



smaller than those obtained through fits of iDraine fe Lil ( 120071 ) models, and the disagreement 



is larger for galaxies with cooler far-infr ared colors. This systematic difference is due to the 



superior ability of the lDraine fe Lil ( 120071 ) dust model to represent the dust emission spectrum 
from the near-infrared through the sub-millimeter, with (for a given value of ^pah) a single 
dust opacity function n u , but allowing for a distribution of starlight heating intensities and 
resulting dust temperature distributions. 



Herschel is an ESA space observatory with science instruments provided by European- 
led Principal Investigator consortia and with important participation from NASA. IRAF, 
the Image Reduction and Analysis Facility, has been developed by the National Optical 
Astronomy Observatories and the Space Telescope Science Institute. 
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Table 1. Galaxy Sample 



Galaxy 


Alternative 


Optical 


E(B — V) 


ao & <5o 


£>25 


2a 


2b 


PA 


TIR 




Name 


Morph. 


(mag.) 


(J2000) 


( ) 


in \ 
( ) 


iii\ 
( ) 


( ) 


IT \ 


NGC0337 




SBd 


0.112 


005950. 7— 073444 


2.9 


253 


194 


140 


10.1 


NGC0584 




E4 


0.042 


01 31 20 6—065205 


4.2 


326 


278 


60 


8.8 


NGO0628 


UGCm 149 


SAc 


0.070 


013642.4+154711 


10.5 
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90 


9.9 


NGC0855 


UGC01718 


E 


0.071 


021403.7+275237 


2.6 


259 


169 


60 


8.6 


NGC0925 




SABd 


0.076 


02271 3 6+333504 


10.5 


735 


486 


105 


9.7 


NGC1097 


UGCA041 


SBb 


0.027 


024618.0—301642 


9.3 


758 


612 


130 


10.7 


NGC1266 




SBO 


0.098 


031600 7—022541 


1.5 


234 


232 





10.4 


NGC1291 




SBO/a 


0.013 


031717.9—410616 


9.8 


884 


836 


90 


9.5 


NGC1316 


FornaxA 


SABO 


0.021 


032241.2—371210 


12.0 


864 


583 


50 


9.9 


NGC1377 




SO 


0.028 


033639.0—205408 


1.8 


181 


162 


90 


10.1 


NGC1404 




El 


0.011 


033852.3—353540 


3.3 


524 


369 


149 




IC0342 


UGC02847 


SABcd 


0.558 


034648.5+680538 


21.4 


1621 


1433 


100 


10.1 


NGC1482 




SAO 


0.040 


035439.0—203009 


2.5 


349 


310 


119 


10.6 


NGC1512 




SBab 


0.011 


040355.6—432149 


8.9 


1001 


928 


83 


9.5 


NGC2146 


TJGO03429 


Sbab 


0.096 


061835 6+782129 


6.0 


236 


235 


120 


11.0 


HoII 




Im 


0.032 


081910 8+704320 


7.9 


554 


465 


60 


7.8 


DDO053 


UGC04459 


Im 


0.038 


083407.4+661043 


1.5 


155 


142 


90 


7.0 


NGC2798 


UGC04905 


SBa 


0.020 


091723 1+415957 


2.6 


235 


232 


90 


10.6 


NGC2841 


UGC04966 


SAb 


0.015 


092203.3+505837 


8.1 


629 


334 


150 


10.1 


NGC2915 




10 


0.275 


092609.4—763736 


1.9 


183 


132 


110 


7.6 


Hoi 


UGC051 39 


IABm 


0.050 


094033 6+71 1 1 20 


3.6 


264 


219 


63 


7.1 


NGC2976 


UGC05221 


SAc 


0.071 


094715.3+675509 


5.9 


541 


353 


144 


8.9 


NGC3049 


UGC05325 


SBab 


0.038 


095449 6+091614 


2.2 


218 


160 


29 


9.5 


NGC3077 


UGC05398 




0.067 


100317 5+684354 


5.4 


488 


436 


64 


8.9 


M81dwB 


UGC05423 


Im 


0.081 


100531.2+702151 


0.9 


134 


90 


139 


6.5 


NGC319D 


UGC05559 


SAap 


0.025 


101805 7+214957 


4.4 


334 


196 


117 


9.9 


NGC3184 


UGC05557 


SABcd 


0.017 


101815.6+412542 


7.4 


614 


538 


169 


10.0 


NGC3198 


UGC05572 


SBc 


0.012 


101954 8+453301 


8.5 


518 


315 


35 


10.0 


IC2574 


UGC05666 


SABm 


0.036 


102823.9+682505 


13.2 


864 


486 


59 


8.3 


NGC3265 


UGC05705 


E 


0.024 


103106.8+284751 


1.3 


184 


175 


50 


9.4 


NGC3351 


M095 


SBb 


0.028 


104358.1+114210 


7.4 


592 


441 


11 


9.9 


NGC3521 


UGC06150 


SABbc 


0.057 


110548.1-000127 


11.0 


926 


455 


165 


10.5 


NGC3621 


UGCA232 


SAd 


0.081 


111818.3-324855 


12.3 


791 


555 


160 


9.9 


NGC3627 


M066 


SABb 


0.033 


112013.4+125927 


9.1 


745 


486 


167 


10.4 


NGC3773 


UGC06605 


SAO 


0.027 


113813.1+120644 


1.2 


118 


116 





8.8 


NGC3938 


UGC06856 


SAc 


0.021 


115250.3+440715 


5.4 


504 


468 





10.3 


NGC4236 


UGC07306 


SBdm 


0.015 


121643.2+692719 


21.9 


1240 


369 


162 


8.7 


NGC4254 


M099 


SAc 


0.039 


121849.7+142519 


5.4 


519 


420 


60 


10.6 


NGC4321 


M100 


SABbc 


0.026 


122254.8+154907 


7.4 


558 


483 


40 


10.5 


NGC4536 


UGC07732 


SABbc 


0.018 


123427.5+021113 


7.6 


454 


376 


120 


10.3 


NGC4559 


UGC07766 


SABcd 


0.018 


123558.1+275752 


10.7 


576 


327 


140 


9.5 


NGC4569 


M090 


SABab 


0.047 


123650.2+131001 


9.5 


593 


327 


21 


9.7 


NGC4579 


M058 


SABb 


0.041 


123743.8+114858 


5.9 


325 


271 


90 


10.1 


NGC4594 


M104 


SAa 


0.051 


123959.6-113726 


8.7 


767 


669 


90 


9.6 


NGC4625 


UGC07861 


SABmp 


0.018 


124154.8+411623 


2.2 


298 


214 


100 


8.8 
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Table 1 — Continued 



Galaxy 


Alternative 


Optical 


E(B - V) 


ao &c So 


£>25 


2a 


2b 


PA 


TIR 




Name 


Morph. 


(mag.) 


(J2000) 


(') 


(") 


(") 


(°) 


(Lq) 


NGC4631 


UGC07865 


SBd 


0.017 


124204.2+323219 


15.5 


901 


240 


85 


10.4 


NGC4725 


UGC07989 


SABab 


0.012 


125027.7+252948 


10.7 


689 


523 


30 


9.9 


NGC4736 


M094 


SAab 


0.018 


125055.2+410652 


11.2 


944 


899 





9.8 


DD0154 


UGC08024 


IBm 


0.009 


125407.6+270916 


3.0 


216 


126 


50 




NGC4826 


M064 


SAab 


0.041 


125643.3+214048 


10.0 


716 


427 


114 


9.6 


DD0165 


UGC08201 


1m 


0.024 


130625.9+674229 


3.5 


263 


161 


90 




NGC5055 


M063 


SAbc 


0.018 


131549.2+420147 


12.6 


1097 


711 


80 


10.3 


NGC5398 


Tololo89 


SBdm 


0.066 


140121.2-330402 


2.8 


198 


146 





8.6 


NGC5457 


M101 


SABcd 


0.009 


140325.0+542429 


28.8 


1800 


1446 


37 


10.4 


NGC5408 




IBm 


0.068 


140321.1-412241 


1.6 


256 


209 


67 


8.3 


NGC5474 


UGC09013 


SAcd 


0.011 


140500.8+533920 


4.8 


412 


373 


90 


8.7 


NGC5713 


UGC09451 


SABbcp 


0.039 


144011.4-001726 


2.8 


225 


225 


90 


10.5 


NGC5866 


UGC09723 


SO 


0.013 


150628.8+554551 


4.7 


500 


306 


129 


9.8 


NGC6946 


UGC11597 


SABcd 


0.342 


203449.2+600959 


11.5 


953 


928 





10.5 


NGC7331 


UGC12113 


SAb 


0.091 


223704.3+342435 


10.5 


683 


335 


168 


10.7 


NGC7793 




SAd 


0.019 


235749.9-323525 


9.3 


716 


526 


98 


9.3 



Note. — D25 is the diameter of the B band isophote defined at 25 mag arcsec -2 . 2a and 2b 
are the lengths of the major and minor axes used in the elliptical aperture photometry described 
herein; the position angle of the aperture's major axis is measured east of nor th. The total infrared 
listed in the last column i s deriv ed from E quation 4 of Dale fc Heloul ([20021 ) and the far- infrared 
fluxes in bale et all |2007l . l2009t ) and from lEngelbracht et al.l (|2008h for IC 0342 and NGC 2146. 
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Table 2. Far-Infrared/Sub-millimeter Flux Densities 



Galaxy 



PACS 
70 

(Jy) 



PACS 
100 /im 

(Jy) 



PACS 
160 (im 

(Jy) 



SPIRE 
250 /urn 

(Jy) 



SPIRE 
350m 

(Jy) 



SPIRE 
500 /im 

(Jy) 



NGC0337 

NGC0584 

NGC0628 

NGC0855 

NGC0925 

NGC1097 

NGC1266 

NGC1291 

NGC1316 

NGC1377 

NGC1404 

IC0342 

NGC1482 

NGC1512 

NGC2146 

HoII 

DDO053 

NGC2798 

NGC2841 

NGC2915 

Hoi 

NGC2976 

NGC3049 

NGC3077 

M81dwB 

NGC3190 

NGC3184 

NGC3198 

IC2574 

NGC3265 

NGC3351 

NGC3521 

NGC3621 

NGC3627 

NGC3773 

NGC3938 

NGC4236 

NGC4254 a 

NGC4321 a 

NGC4536 a 

NGC4559 

NGC4569 a 

NGC4579 a 

NGC4594 



1.30±0.07E+1 1.95±0.10E+1 1.96±0.10E+1 



3.67±0.18E+1 
2.30±0.12E+0 
1.08±0.06E+1 
7.75±0.39E+1 
1.45±0.07E+1 
5.26±0.32E+0 
5.81±0.33E+0 
6.89±0.35E+0 
<6.71E-1 
4.48±0.22E+2 
4.07±0.20E+1 
7.99±0.47E+0 
1.98±0.10E+2 
3.18±0.35E+0 

3.90±0.42E-1 
2.42±0.12E+1 
9.49±0.49E+0 
1.01±0.06E+0 

3.71±0.62E-1 
1.92±0.10E+1 
3.40±0.18E+0 
2.04±0.10E+1 

1.21±0.41E-1 
6.30±0.33E+0 
1.55±0.08E+1 
9.75±0.51E+0 
5.61±0.37E+0 
2.47±0.13E+0 
2.53±0.13E+1 
7.85±0.39E+1 
4.95±0.25E+1 
1.04±0.05E+2 
1.29±0.08E+0 
1.58±0.08E+1 
7.46±0.46E+0 
5.64±0.28E+1 
4.12±0.21E+1 
3.89±0.20E+1 
1.59±0.08E+1 
1.46±0.07E+1 
9.94±0.51E+0 
7.87±0.49E+0 



7.40±0.37E+1 
2.04±0.12E+0 
2.47±0.12E+1 
1.16±0.06E+2 
1.59±0.08E+1 
1.28±0.07E+1 
9.30±0.50E+0 
5.85±0.30E+0 
<7.06E-1 
8.47±0.42E+2 
4.95±0.25E+1 
1.38±0.07E+1 
2.32±0.12E+2 
3.89±0.45E+0 

4.80±1.21E-1 
2.73±0.14E+1 
2.57±0.13E+1 
1.66±0.09E+0 

4.21±0.70E-1 
3.58±0.18E+1 
4.59±0.23E+0 
2.79±0.14E+1 

2.01±0.31E-1 
1.06±0.05E+1 
3.47±0.17E+1 
2.00±0.10E+1 
7.57±0.42E+0 
3.10±0.16E+0 
4.61±0.23E+1 
1.58±0.08E+2 
9.44±0.47E+1 
1.79±0.09E+2 
1.85±0.11E+0 
2.86±0.15E+1 
1.23±0.07E+1 
1.06±0.05E+2 
8.55±0.43E+1 
5.26±0.26E+1 
3.10±0.16E+1 
3.04±0.15E+1 
2.33±0.12E+1 
2.39±0.13E+1 



1.16±0.06E+2 
2.16±0.12E+0 
3.65±0.18E+1 
1.34±0.07E+2 
1.13±0.06E+1 
2.03±0.11E+1 
1.15±0.06E+1 
3.38±0.19E+0 
<1.02E 
l.ll±0.06E+3 
4.20±0.21E+1 
1.87±0.10E+1 
1.81±0.09E+2 
3.86±0.63E+0 

2.50±1.77E-1 
2.06±0.10E+1 
4.95±0.25E+1 
1.46±0.11E+0 

3.72±1.20E-1 
4.64±0.23E+1 
4.54±0.24E+0 
2.83±0.14E+1 

2.42±0.82E-1 
1.54±0.08E+1 
5.49±0.28E+1 
2.99±0.15E+1 
9.61±0.53E+0 
2.63±0.15E+0 
5.51±0.28E+1 
2.10±0.10E+2 
1.28±0.06E+2 
2.02±0.10E+2 
1.91±0.14E+0 
3.96±0.20E+1 
1.85±0.11E+1 
1.30±0.07E+2 
1.20±0.06E+2 
5.55±0.28E+1 
4.10±0.21E+1 
4.03±0.20E+1 
3.54±0.18E+1 
3.89±0.20E+1 



9.79±0.70E+0 
<8.84E-1 
6.55±0.47E+1 
1.48±0.11E+0 
2.77±0.20E+1 
7.22±0.51E+1 
4.38±0.31E+0 
1.59±0.11E+1 
4.80±0.37E+0 
1.32±0.10E+0 
<4.30E-1 
5.95±0.42E+2 
1.68±0.12E+1 
1.56±0.11E+1 
6.55±0.47E+1 
1.82±0.16E+0 

1.86±0.32E-1 
8.02±0.57E+0 
3.49±0.25E+1 

9.28±0.73E-1 

3.56±0.53E-1 
2.50±0.18E+1 
2.80±0.20E+0 
1.43±0.10E+1 

1.87±0.25E-1 
8.88±0.63E+0 
3.43±0.24E+1 
1.96±0.14E+1 
7.16±0.52E+0 
1.24±0.10E+0 
3.24±0.23E+1 
1.14±0.08E+2 
7.12±0.51E+1 
9.67±0.69E+1 
1.06±0.08E+0 
2.37±0.17E+1 
1.16±0.08E+1 
6.57±0.47E+1 
6.76±0.48E+1 
2.88±0.20E+1 
2.55±0.18E+1 
2.24±0.16E+1 
2.12±0.15E+1 
2.56±0.18E+1 



4.37±0.31E+0 
<8.19E-1 
3.06±0.22E+1 

7.47±0.65E-1 
1.48±0.11E+1 
3.08±0.22E+1 
1.60±0.12E+0 
7.98±0.59E+0 
2.06±0.19E+0 

4.94±0.47E-1 
<3.99E-1 
2.61±0.19E+2 
6.35±0.45E+0 
8.66±0.64E+0 
2.33±0.17E+1 

8.04±1.08E-1 

9.99±2.81E-2 
2.90±0.21E+0 
1.60±0.11E+1 

5.28±0.47E-1 

2.23±0.47E-1 
1.17±0.08E+1 
1.41±0.11E+0 
6.47±0.47E+0 

1.03±0.22E-1 
3.71±0.27E+0 
1.53±0.11E+1 
9.95±0.71E+0 
4.83±0.36E+0 

5.51±0.51E-1 
1.37±0.10E+1 
4.72±0.34E+1 
3.17±0.23E+1 
3.76±0.27E+1 

4.34±0.38E-1 
1.03±0.07E+1 
7.37±0.54E+0 
2.66±0.19E+1 
2.79±0.20E+1 
1.26±0.09E+1 
1.28±0.09E+1 
9.41±0.67E+0 
8.87±0.63E+0 
1.21±0.09E+1 



1.88±0.14E+0 
<7.63E-1 
1.33±0.10E+1 

2.66±0.40E-1 
8.03±0.58E+0 
1.26±0.09E+1 

5.32±0.55E-1 
3.52±0.29E+0 

8.16±1.33E-1 

1.99±0.32E-1 
<3.72E-1 
1.02±0.07E+2 
2.21±0.17E+0 
4.20±0.34E+0 
7.45±0.53E+0 

3.37±1.69E-1 
<1.25E-1 
1.08±0.09E+0 
7.01±0.50E+0 

2.54±0.32E-1 

1.35±0.41E-1 
4.79±0.35E+0 

7.97±0.65E-1 
2.89±0.22E+0 

5.66±2.83E-2 
1.38±0.11E+0 
6.73±0.49E+0 
4.74±0.34E+0 
2.13±0.19E+0 

2.38±0.35E-1 
5.32±0.39E+0 
1.94±0.14E+1 
1.46±0.10E+1 
1.44±0.10E+1 

1.80±0.24E-1 
4.34±0.32E+0 
4.21±0.32E+0 
9.16±0.66E+0 
1.08±0.08E+1 
5.53±0.40E+0 
6.37±0.46E+0 
3.67±0.27E+0 
3.54±0.26E+0 
5.56±0.41E+0 
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Table 2 — Continued 



Galaxy 




PACS 




PACS 




PACS 


SPIRE 


SPIRE 


SPIRE 






70 /im 




100 /im 




160 (im 


250 fj,m 


350m 


500 /im 






(Jy) 




(Jy) 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


NGC4625 


1 


36±0.12E+0 


3 


04±0.20E+0 


1 


48±0.23E+0 


2.81±0.21E+0 


1.40±0.11E+0 


6.44±0.62E-1 


NGC4631 


1 


37±0.07E+2 


2 


23±0.11E+2 


2 


46±0.12E+2 


1.24±0.09E+2 


5.45±0.39E+1 


2.40±0.17E+1 


NGC4725 a 


7 


93±0.46E+0 


2 


28±0.12E+1 


4 


66±0.23E+1 


3.27±0.23E+1 


1.66±0.12E+1 


7.93±0.57E+0 


NGC4736 


1 


03±0.05E+2 


1 


59±0.08E+2 


1 


45±0.07E+2 


7.04±0.50E+1 


2.80±0.20E+1 


1.18±0.09E+1 


DD0154 




<3.31E-1 




<4.27E-1 




<2.67E-1 


<1.62E-1 


<1.50E-1 


<1.39E-1 


NGC4826 


5 


47±0.27E+1 


9 


57±0.48E+1 


9 


41±0.47E+1 


4.24±0.30E+1 


1.64±0.12E+1 


6.30±0.46E+0 


DD0165 




<4.13E-1 




<5.33E-1 




<3.33E-1 


<2.01E-1 


<1.87E-1 


<1.74E-1 


NGC5055 


7 


34±0.37E+1 


1 


70±0.08E+2 


2 


48±0.12E+2 


1.50±0.11E+2 


6.42±0.46E+1 


2.61±0.19E+1 


NGC5398 


2 


19±0.12E+0 


2 


98±0.16E+0 


2 


75±0.16E+0 


2.03±0.15E+0 


1.05±0.08E+0 


5.52±0.49E-1 


NGC5457 


1 


23±0.06E+2 


2 


43±0.12E+2 


3 


41±0.17E+2 


2.08±0.15E+2 


9.69±0.69E+1 


4.53±0.32E+1 


NGC5408 


3 


60±0.19E+0 


2 


65±0.15E+0 


2 


02±0.11E+0 


7.85±0.72E-1 


3.86±0.50E-1 


2.09±0.42E-1 


NGC5474 


3 


24±0.18E+0 


1 


61±0.25E+0 


7 


12±0.37E+0 


5.37±0.39E+0 


2.91±0.22E+0 


1.58±0.13E+0 


NGC5713 


2 


89±0.14E+1 


4 


03±0.20E+1 


3 


93±0.20E+1 


1.68±0.12E+1 


6.39±0.46E+0 


2.30±0.17E+0 


NGC5866 


8 


12±0.42E+0 


1 


67±0.09E+1 


1 


84±0.10E+1 


8.04±0.58E+0 


3.14±0.23E+0 


1.14±0.10E+0 


NGC6946 


2 


46±0.12E+2 


1 


35±0.22E+2 


5 


42±0.27E+2 


2.74±0.19E+2 


1.09±0.08E+2 


4.28±0.30E+1 


NGC7331 


6 


53±0.33E+1 


1 


32±0.07E+2 


1 


76±0.09E+2 


9.53±0.68E+1 


4.06±0.29E+1 


1.65±0.12E+1 


NGC7793 


3 


20±0.16E+1 


6 


58±0.33E+1 


9 


11±0.46E+1 


5.63±0.40E+1 


2.84±0.20E+1 


1.39±0.10E+1 



Note. — The compact table entry format T.UViW.XYEZ implies (T.UV±W.XY)x 10 z . See §[3] 
for corrections that have been applied to the data. The uncertainties include both statistical and 
systematic effects. 5er upper limits are provided for non-detections. No color corrections have been 
applied. 

a SPIRE imaging taken from the Herschel Reference Survey |Boselli et al.lboid) . 
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Fig. 1 — A PACS 160 /an image of KINGFISH galaxy NGC 7331. The large red ellipse 
indicates the photometric aperture listed in Table [fl the green circles show the sky apertures, 
and the two blue circles identify sources to be removed before the photometry is executed. 
North is up, East is to the left. 
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Fig. 2. — A comparison of spatially-integrated Spitzer/MIPS and Herschel/PACS photome- 
try at 70 /im (top) and 160 /zm (bottom) for all sample galaxies. The horizontal dotted lines 
indicate the ratios corresponding to perfect agreement between data taken by the two obser- 
vatories, after accounting for differences in the Spitzer and Herschel spectral responses and 
calibration schemes; for typical galaxy spectral energy distributions, the dotted lines differ 
from unity by a factor of 1.06 and 1.015 for the 70 and 160 /im comparisons, respectively. 
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Fig. 3. — A far- infrared/sub-millimeter color-color diagram for the KINGFISH sam ple. 
"Chara cteristic" oxygen abundan ces are taken from Table 9 of iMoustakas et al.l (120101 ). us- 
ing the iPilyugin &: Thuanl (120051 ) metallicity scale; if u navailable in Tabl e 9, th en the value 
is computed using the B luminosity and Equation 10 of IMoustakas et al.l (120101 ). which uses 
the same metal licity scale. The soli d line indicates the sequence of model spectral energy 
distributions of iDale fc Heloul (120021 ). derived from the average global trends for a sample 
of normal star-forming galaxies observed by ISO and IRAS. A set of typical error bars is 
provided. 
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Fig. 4. — A second far- infrared/ sub-millimeter color-color diagram for the KINGFISH sample 
(see also Figure [3]). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for all the 
galaxies in the KINGFISH sample, sorted by Right Ascension. Herschel data are represented 
by filled circles and ancillary data are indicated by open symbols (triangles: 2MASS and 
IRAS; circles: Spitzer, squares: ISO and SCUBA). Arrows indicate ha upper limits. The 
solid curve is the sum of a 5000 K stellar blackbody (short dashed) along with models of dust 
emission from PDRs (dotted; U > U m \ n ) an d the diffuse i nterst ellar medium (long dashed; 
U = U m i n ). The fitted parameters from these iDraine fc Lil (120071 ) model fits are listed within 
each panel along with the reduced \ 2 ( see § 14.41 for details). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 5. — Globally-integrated infrared/sub- millimeter spectral energy distributions for the 
KINGFISH sample (continued). 
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Fig. 6. — The collection of iDraine &: Lil (12007) model fits, color c oded according to star 
formation rate in M & yr -1 in the top panel (jKennicutt et al.l 120121 1 and log 10 (LTm/ X ) in 
the bottom panel (see Tabled]). 
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Fig. 7. — Comparison of the dust model parameters obtained from fitting iDraine fc Lil (120071 ) 
spectral energy distribution models to the entire observed photometric data from 3.6 to 
500 /iin, for fits including both Spitzer and Herschel data ('S+H') versus fits using just the 
Spitzer photometry ('S'). The parameters constrain the quantity of interstellar dust and 
their heating; see § 14.31 The comparison is made as a function of the / I/ (70//m) / f v (160fim) 
ratio, which is related to the average temperature of the interstellar dust grains. Reference 
dotted lines are drawn for a ratio of unity. 
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Fi g. 8. — The same a s Figu re [7] except plotted as a func tion of oxygen abundance as presented 
in iMoustakas et al.l (120101 ) . using data placed on the iPilyugin fc Thuan metallicity 
scale. 
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Fig. 9. — The ratio of single-temperature, blackbody-based dust mass to that obtained 
from lDraine fc Lil (120071 ) model fits to the observed infrared/sub-millimeter spectral energy 
distributions. Top (Bottom): The modified blackbody fits are based on the (error- weighted) 
Herschel 100, 160, 250, 350, and 500 /im (70, 100, 160, 250, 350, and 500 fim) photometry, 
and both Td and the dust opacity coefficient (3 are allowed to freely vary. Reference dotted 
lines are drawn for a ratio of unity. 
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Fig. 10. — The distributions of dust temperature and dust emissivity index (5 when single- 
temperature modified blackbodies are fit to the (100-500 /im) infrared spectral energy dis- 
tributions. 
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Fig. 11. — A comparison of fitting iDraine fc Lil (120071 ) models to 3.6-500 /im data versus fits 
of single temperature blackbodies to 100-500 /im photometry, for a galaxy with warm dust 
(NGC 2798) and a galaxy with cool dust (NGC 2841). 



